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Abstract: On directed Small-World networks the Majority-vote model with 
noise is now studied through Monte Carlo simulations. In this model, the order- 
disorder phase transition of the order parameter is well defined in this system. We 
calculate the value of the critical noise parameter q c for several values of rewiring 
probability p of the directed Small- World network. The critical exponentes (3/v, 
^ jv and 1/u were calculated for several values of p. 
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Introduction 

Sumour and Shabat pp, [2] investigated Ising models on directed Barabasi- 
Albert networks [3] with the usual Glauber dynamics. No spontaneous mag- 
netisation was found, in contrast to the case of undirected Barabasi- Albert 
networks [U El E] where a spontaneous magnetisation was found lower a 
critical temperature which increases logarithmically with system size. Lima 
and Stauffer [7] simulated directed square, cubic and hypercubic lattices in 
two to five dimensions with heat bath dynamics in order to separate the 
network effects form the effects of directedness. They also compared differ- 
ent spin flip algorithms, including cluster flips [SJ, for Ising-Barabasi-Albert 
networks. They found a freezing-in of the magnetisation similar to [H [2], 
following an Arrhenius law at least in low dimensions. This lack of a spon- 
taneous magnetisation (in the usual sense) is consistent with the fact that if 
on a directed lattice a spin Sj influences spin Si, then spin Si in turn does 
not influence Sj, and there may be no well-defined total energy. Thus, they 
show that for the same scale-free networks, different algorithms give different 
results. It has been argued that nonequilibrium stochastic spin systems on 
regular square lattice with up-down symmetry fall in the universality class of 
the equilibrium Ising model [9]. This conjecture was found in several models 
that do not obey detailed balance pH [TE1 H21 H31 CH] • Lima [151 CE] inves- 
tigated the majority-vote model on directed and undirected Barabasi- Albert 
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network and calculated the ft/v, j/u, and \ jv exponents and these are dif- 
ferent from the Ising model and depend on the values of connectivity z of the 
Barabasi-Albert network. Campos et al. [T7] investigated the majority-vote 
model on undirected small-world network by rewiring the two-dimensional 
square lattice. These small-world networks, aside from presenting quenched 
disorder, also possess long-range interactions. They found that the critical 
exponents 7/V and (3/v are different from the Ising model and depend on 
the rewiring probability. However, it was not evident whether the exponent 
change was due to the disordered nature of the network or due to the pres- 
ence of long-range interactions. Lima et al. [18j studied the majority-vote 
model on Voronoi-Delaunay random lattices with periodic boundary condi- 
tions. These lattices posses natural quenched disorder in their conecctions. 
They showed that presence of quenched connectivity disorder is enough to 
alter the exponents (3/v and 7/V the pure model and therefore that is a 
relevant term to such non-equilibrium phase-transition. Now, we calculate 
the same (3/v, 7/V, and 1/V exponents for majority- vote model on directed 
small- world networks of Sanchez et al. [T9] , and for these networks the ex- 
ponents are different from the two-dimensional Ising model and independent 
on the values of rewiring probability p of the directed small-world networks. 
Here we study the cases for p > 0, but we verify that when p = our results 
are agree with results of J. M. Oliveira [T3] . 



Model and Simulaton 

We consider the majority-vote model, on directed small-world networks, 
defined [T3J, EDI ESI EI] by a set of "voters" or spins variables a taking the 
values +1 or —1, situated on every site of an directed small- world networks 
with N = L x L sites, were L is the side of square lattice, and evolving in 
time by single spin-flip like dynamics with a probability W{ given by 



where S(x) is the sign ±1 of x if x 7^ 0, S(x) = if x = 0, and the sum runs 
over four neighbours of Cj. In this network, created for Sanchez et al. [19], 
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Figure 1: Sketch of a directed small-world networks constructed from a square 
regular lattice in d — 2. Figure gentily yielded by Juan M. Lopez of the paper 
of Sanchez et al. [19] . 

see Fig. 1, we start from a two-dimensional square lattice consisting of sites 
linked to their four nearest neighbors by both outgoing and incoming links. 
Then, with probability p, we reconnect nearest-neighbor outgoing links to a 
different site chosen at random. After repeating this process for every link, we 
are left with a network with a density p of SW directed links. Therefore, with 
this procedure every site will have exactly four outgoing links and a varying 
(random) number of incoming links. The control parameter q plays the role 
of the temperature in equilibrium systems and measures the probability of 
aligning antiparallel to the majority of neighbours. 

To study the critical behavior of the model we define the variable m = 
Sili cri/N. In particular, we were interested in the magnetisation, suscepti- 
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bility and the reduced fourth-order cumulant: 



M(q) = [< \m\ >} av , 
X(q) = N[< m 2 > - < \m\ > 2 } 



(2) 



(3) 



17(g) = 1 



< m 4 > 



(4) 



3 < \m\ > 2 . 



av 



where < ... > stands for a thermodynamics average and [...] a „ square brackets 
for a averages over the 20 realizations. 

These quantities are functions of the noise parameter q and obey the 
finite-size scaling relations 



where u, (3, and 7 are the usual critical exponents, fi(x) are the finite size 
scaling functions with 



being the scaling variable, and the brackets [1 + ...] indicate corretions-to- 
scaling terms. Therefore, from the size dependence of M and x we obtained 
the exponents (3/v and 7///, respectively. The maximum value of susceptibil- 



were the constant b is close to unity. Therefore, the relations (7) and (9) are 
used to determine the exponent \jv. 

We have performed Monte Carlo simulation on directed SW network with 
various values of probability p. For a given p, we used systems of size L — 8, 
16, 32, 64, and 128. We waited 10000 Monte Carlo steps (MCS) to make the 
system reach the steady state, and the time averages were estimated from 
the next 10000 MCS. In our simulations, one MCS is accomplished after all 
the iV spins are updated. For all sets of parameters, we have generated 20 
distinct networks, and have simulated 20 independent runs for each distinct 
network. 



M = L-V v f m (x)[l + ...], 
X = L^f x (x)[l + ...], 

%- = L 1 "'f u (x)[l + ...], 



(5) 
(6) 



(7) 




(8) 
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Results and Discussion 

In Fig. 2 we show the dependence of the magnetisation M and the sus- 
ceptiblity \ on the noise parameter, obtained from simulations on directed 
SW network with L = 128 x 128 sites and several values of probability p. In 
part (a) each curve for M, for a given value of L and p, suggests that there 
is a phase transition from an ordered state to a disordered state. The phase 
transition occurs at a value of the critical noise parameter q c , which is an 
increasing function the probability p of the directed SW network. In part (b) 
we show the corresponding behavior of the susceptibility x> the value of q 
where \ has a maximum is here identified as q c . In Fig. 3 we plot Binder's 
fourth-order cumulant for different values of L and two different values of p. 
The critical noise parameter q c , for a given value of p, is estimated as the 
point where the curves for different system sizes L intercept each other. In 
Fig 4 the phase diagram is shown as a function of the critical noise parameter 
q c on probability p obtained from the data of Fig. 3. 

The phase diagram of the majority-vote model on directed SW network 
shows that for a given network (fixed p ) the system becomes ordered for 
q < q c , whereas it has zero magnetisation for q > q c . We notice that the 
increase of q c as a function the p is slower that the one in |21j. In Figs. 5 and 
6 we plot the dependence of the magnetisation and susceptibility, respectively, 
at q = q c versus the system size L. The slopes of curves correspond to the 
exponent ratio (3/v and 7/1/ of according to Eq. (5) and (6), respectively. 
The results show that the exponent ratio (3/v and 7/1/ at q c are independent 
of p (along with errors), see Table I. 

In Fig. 7 we display the scalings for susceptibility at q — q c {L) (square), 
x{q c {L)), and for its maximum amplitude, Xi] ax ', and the scalings for suscep- 
tibility at the q = q c obtained from Binder's cumulant, x(Qc) (circle), versus 
L for probability p = 0.5. The exponents ratio 7/V are obtained from the 
slopes of the straight lines. For almost all the values of p, the exponents 
7/// of the two estimates agree (along with errors). We also observe that an 
increased p does not mean a tendency to increase or decrease the exponent 
ratio 7/V, see Table I. Therefore we can use the Eq. (9) , for fixed p, to 
obtain the critical exponent 1/v, see Fig. 8. 

To improve our results obtained above we start with all spins up, a number 
of spins equal to N = 640000, and time up 2,000,000 (in units of Monte 
Carlo steps per spins). Then we vary the noise parameter q and at each q 
study the time dependence for 9 samples. We determine the time r after 
which the magnetisation has flipped its sign for fisrt time, and then take the 
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median values of our nine samples. So we get different values t% for different 
noise parameters q. In Fig. 9 show that the decay time goes to infinity at 
some q positive values this behavior sure that there is a phase transition for 
Majority- vote on directed SW network. 

The Table I summarizes the values of q c , the exponents /3/u, and \jv 
. J. M. Oliveira [13] showed that the majority- vote model defined on a regular 
lattice has critical exponents that fall into the same class of universality as 
the corresponding equilibrium Ising model. Campos et al [17] investigated 
the critical behavior of the majority- vote on small-world networks by rewiring 
the two-dimensional square lattice, Pereira et al [21] studied this model on 
Erdos-Renyi's random graphs, and Lima et al [TB] also studied this model 
on Voronoy-Delaunay lattice and Lima on directed Barabasi- Albert network 
[T5] . The results obtained these authors show that the critical exponents of 
majority-vote model belong to different universality classes. 



V 


Qc 


P/u 


7 /zA 




l/i/ 


0.1 


0.122(3) 


0.423(17) 


1.178(13) 


1.214(39) 


0.837(223) 


0.3 


0.149(3) 


0.419(21) 


1.148(5) 


1.152(28) 


1.059(208) 


0.5 


0.160(2) 


0.441(12) 


1.116(5) 


1.120(25) 


1.010(52) 


0.8 


0.164(2) 


0.436(9) 


1.149(5) 


1.117(23) 


1.248(158) 


1.0 


0.165(2) 


0.415(18) 


1.139(8) 


1.122(25) 


1.032(81) 



Table 1: The critical noise q c , and the critical exponents , for directed SW 
network with probability p. Error bars are statistical only. 

Finally, we remark that our MC results obtained on directed SW network 
for majority-vote model show that critical exponents are different from the 
results of [13] for regular lattice, of Pereira et al [21] for Erdos-Renyi's random 
graphs, Lima [15] and Campos et al. [17] for on undirected SW network. 

Conclusion 

In conclusion, we have presented a very simple nonequilibrium model 
on directed SW network [19]. In these networks, the majority- vote model 
presents a second-order phase transition which occurs with probability p > 0. 
The exponents obtained are different from the other models [T3J HH H51 EHl 
[UJET] suggesting that these exponents belong to another class of universality. 
However, the exponents in the critical point q c , [3/u, 7/1/ and \jv when p 
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grows no increase and do not decrease and are independents of a growing p 
for p > 0. 
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Magnetisation versus q for p=0.1, 0.3, 0.5, 0.8 &1.0 
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Figure 2: Magnetisation and susceptibility as a function of the noise param- 
eter q, for N = 16384 sites. From left to right we have p — 0.1, 0.3, 0.5, 0.8, 
and 1.0 . 
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Figure 3: Binder's fourth-order cumulant as a function of q. We have p = 0.3 
and p = 0.8 for L = 8, 16, 32, 64 and 128. 
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Figure 4: The phase diagram, showing the dependence of critical noise 
rameter q c on probability p. 
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Figure 6: In x(Qc) versus In L. From top to bottom p = 0.1, 0.3, 0.5, 0.8, 
and 1.0. 
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directed SW network: p=0.5 and 640,000 sites 
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Figure 9: Plot of l/ln(tau) versus q for p=0.5. 
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